Phase diagram for quantum Hall states in graphene 
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We investigate integer and half-integer filling states (uniform and unidimensional stripe states 
respectively) for graphene using the Ifartree-Fock approximation. For fixed filling factor, the ratio 
between the scales of the Coulomb interaction and Landau level spacing g = {e^ / et) / [hvp / 1), with 
I the magnetic length, is a field-independent constant. However, when B decreases, the number 
of filled negative Landau levels increases, which surprisingly turns out to decrease the amount of 
Landau level mixing. The resulting states at fixed filling factor u (for v not too big) have very 
little Landau level mixing even at arbitrarily weak magnetic fields. Thus in the density-field phase 
diagram, many different phases may persist down to the origin, in contrast to the more standard two 
dimensional electron gas, in which the origin is surrounded by Wigner crystal states. We demonstrate 
that the stripe amplitudes scale roughly as _B, so that the density waves "evaporate" continuously 
as i? — > 0. Tight-binding calculations give the same scaling for stripe amplitude and demonstrate 
that the effect is not an artifact of the cutoff procedure used in the continuum calculations. 

PACS numbers: 73.20.Qt, 73.43.-f, 81.05.Uw 



I. INTRODUCTION 

Graphene, a two-dimensional honeycomb lattice of car- 
bon atoms, has attracted intense attention in the past 
few years. ^ Its properties bear some similarities with, 
and some striking differences from, conventional 2D elec- 
tron gas (2DEG) systems found in semiconductor het- 
erostructures. It is well-known that the latter has a rich 
phase diagram in the quantum Hall regime. When r^, 
the average inter-electron distance measured in units of 
Bohr radius, is not very big, there are integer and frac- 
tional quantum Hall liquid states, as well as charge den- 
sity waves (CDWs) of various forms, including Wigner 
crystals of quasi-electrons, bubbles and stripes at fillings 
around these liquid states, and analogous particle-hole 
conjugates of these states^i'^i'^. In high magnetic fields, 
the particular state is essentially determined by the fill- 
ing factor V, defined as the ratio of the electron density 
to the density of magnetic flux quanta penetrating the 
plane. When is increased, these quantum Hall phases 
undergo transitions to Wigner crystal (WC) states with 
a single electron per unit cell. (For very small v, there 
may also be Wigner crystals of composite fermions^i^.) 
If the phase diagram is plotted in the n (density) - B 
(magnetic field) plane, away from the origin, there is a 
fan of quantum Hall phases, but the origin is expected to 
be completely surrounded by Wigner crystal states^ [see 

Fig. [1(5] • 

The integer quantized Hall effect has been observed in 
graphene^iSiiiSiii, and, except for a well-understood shift 
in the precise values of the plateaus^ii^, the Hall conduc- 
tance appears rather similar to that found in the conven- 
tional 2DEG. Nevertheless, the behavior of clean and cold 
graphene in the low doping limit is likely to be different 
than that of the conventional 2DEG. Unlike the latter, 
non-interacting electrons in graphene to a good approx- 



imation obey a massless Dirac equationi^iilii^, with two 
inequivalent Dirac points in two different valleys (denoted 
as K and K') in the Brillouin zone. When the system is 
undoped the Fermi energy passes directly through these 
Dirac points. With interactions, continuum^® and tight- 
bindingii studies of the this system in mean-field theory 
indicate that the system remains in a liquid state in zero 
magnetic field even at arbitrarily low doping. On the 
other hand, Hartree-Fock calculationsi^ and exact diag- 
onalization studies^ suggest that CDWs are possible in 
a large magnetic field - where states are restricted to a 
single or two2S Landau levels (LLs) - and that the phase 
diagram is similar to that of the conventional 2DEG. In 
this paper we address the question of how the system 
passes from these strong-field states into the liquid state 
as the field and density are lowered to small values. 

For the conventional 2DEG, the quantum Hall states 
give way to the WC in the low-field, low density limit 
due to Landau level mixing (LLM). This allows the elec- 
trons to form wavepackets that are more localized than 
is possible within a single Landau level, thereby lowering 
the interaction energy^^. The degree of LLM is deter- 
mined by a coupling constant g, the ratio of the typi- 
cal Coulomb interaction energy Ec to the scale of the 
LL separation. For both graphene and the conventional 
2DEG, Ec is given by e^/{et), where I = ^hc/eB is 
the magnetic length. However, the LL separations are 
different in the two cases. In the conventional 2DEG, it 
is given by hjjc — heB/mc, so that g cx and in 

the large B limit where g is small, LLM is negligible. In 
graphene, the LLs are not equally spaced^^, so we instead 
characterize it by the gap between the n = and n = 1 
LLs divided by ^/2, hvF/£, where vp is the Fermi velocity. 
Then g = {e^/e)/{hvF) is a field-independent constant 
typically estimated to be of order 1 or smaller. Never- 
theless, even though g is field- independent, the degree 
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FIG. 1: Schematic phase diagrams for conventional 2DEG and graphene in the integer quantum Hall regime. Here n is the 
electron density. "2D crystals" referred to in diagrams include bubble, quasiparticle and quasihole states, whose lattice constants 
are determined by the filling factor and magnetic length. These differ qualitatively from the Wigner crystal state where the 
lattice constant is set by the electron density. Fractional quantum Hall states, expected to appear at low filling factors in both 
diagrams, are not shown. The major difference between the two phase diagrams is that for the conventional 2DEG, the origin 
is completely surrounded by the Wigner crystal state, while in the graphene case, many different phases persist down to the 
origin. 



of LLM can change with B even for fixed filling factor, 
and we shall see below that it in fact does, albeit by a 
small amount. This is because in addition to positive 
energy levels, the Dirac equation in a magnetic field sup- 
ports negative energy Landau level states, as well as a 
zero energy LL, for each spin, provided the Zeeman en- 
ergy is neglected. Moreover, the low energy theory of 
graphene involves two such Dirac points (K and K' val- 
leys), so there are two copies of these energy levels in 
the spectrum. When undoped, all the negative energy 
states are filled, as well as one of the two zero-energy 
states^. Added electrons interact with the electrons in 
these filled levels, which changes the effective energy of 
the higher LLs. Because the Landau level structure of 
these filled levels varies with field, the splitting between 
the n — and n — 1 energy levels does not precisely the 
follow the ^/B behavior discussed above. 



In a continuum description, the filling of the negative 
levels is characterized by a (negative) integer ric, which 
denotes the lowest LL which must be filled to accommo- 
date one electron per atom, the density of mobile elec- 
trons of undoped graphene^^'^^. An extra field depen- 



dence thus enters the problem through Uc, and is given 

by 



25/(V3aV2) 
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(1) 



where S is the area of the sample, a = 0.246 nm is the 
lattice constant of the triangular (Bravais) lattice and 
the factor of 4 in the denominator comes from the spin 
and valley degeneracies. The interaction of electrons with 
those in the negative levels is most easily described in the 
Hartree-Fock (HF) approximation, where it appears as a 
contribution to the exchange self-energy. For uniform 
liquid states, we find that the Coulomb energy decreases 
faster with decreasing B than the difference in effective 
energy between the highest occupied level and the lowest 
unoccupied level, so that the ratio between kinetic and 
potential energy actually increases with decreasing B, 
because jnd increases. We will demonstrate a similar 
effect for stripe states, and believe it should be ubiquitous 
for charge-ordered and liquid quantum Hall states. 

Because this effect is a result of filling jrid negative 
LLs, it is a concern that it may be an artifact of the 
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cutoff procedure used in our HF calculations. To check 
this, we performed an analogous calculation for interact- 
ing electrons in a tight-binding model, where no artificial 
cutoff needs to be introduced. We obtain results from 
this model that are very similar to those of the contin- 
uum model. 

The consequence of this is that, for states where LLM 
is small at large values of -B, we expect it remain small, 
and even decrease, with decreasing B. While it is not 
immediately obvious that with g ~ 1 one should find 
weak LLM in these quantum Hall states, this does ap- 
pear to be the case for WC and bubble states^, and 
as we demonstrate below, for stripe and uniform liquid 
states. The surprising result is that, within the Hartree- 
Fock approach, one expects these states to persist to ar- 
bitrarily small field. Thus, many different states persist 
down to the origin of the phase diagram in n - B plane 
[see Fig. l(b)| . Because these states follow trajectories 
of fixed ly in the plane, the density of electrons partic- 
ipating these CDW states decreases with B, such that 
their amplitude scales roughly as B, and the wavelength 
as 1/^/B. The stripe states, and by analogy other CDW 
states, disappear continuously as i? ^ 0, eventually be- 
coming indistinguishable from a uniform liquid state in 
the low field limit. Nevertheless, in principle, for a clean, 
undeformed graphene system, this implies that in prin- 
ciple many states emanate from the _B = n = point in 
the phase diagram. 

We note that this behavior is very specific to the l/r 
form for Coulomb interactions that is natural in this sys- 
tem. For shorter range interactions, where a length scale 
other than the magnetic length is involved in the inter- 
action range, the effective value of g will increase with 
decreasing field as in the standard 2DEG, at low den- 
sities and fields LLM should destabilize the high field 
states, and a WC state should result. Such a situation 
could arise if a metallic gate is sufficiently close to the 
graphene plane to effectively screen the long-range com- 
ponent of the Coulomb interaction. 

More generally, the behavior discussed here may be un- 
derstood as being a consequence of the marginal nature 
of the l/r Coulomb interaction in undoped graphene. As 
has been shown by us elsewhere^, the energy difference 
between Landau levels near the Fermi energy is increased 
by the filled Fermi sea, by an amount proportional to 
log |nc|. This logarithmic increase of the LL spacing with 
increasing jnd can be reinterpreted as the Fermi velocity 
being renormalized upwards as the high-energy cutoff of 
the theory is increased?^. That the LL spacing increases 
slightly with increasing \nc\ as the doping is decreased is 
consistent with interactions being marginally irrelevant 
in this system^. Had it decreased instead, the inter- 
actions would be marginally relevant, one would expect 
to find a WC state near the origin of the n — B phase 
diagram. 

This paper is organized as follows. In Section|TTl we de- 
scribe the continuum limit Hamiltonian, and the Hartree- 
Fock approximation used to study quantum Hall states 



in the presence of LLM. In Section IIIIl we discuss the 
results of these continuum calculations. In Section IIV[ 
we introduce a tight-binding model with Hubbard inter- 
actions, and demonstrate that the suppression of LLM 
found in the continuum calculations is not an artifact of 
our cutoff procedure. Finally, we conclude with a sum- 
mary in Section |Vl 



II. HARTREE-FOCK FOR CONTINUUM 
MODEL 

In standard 2DEG's, it is known that the Hartree-Fock 
approximation is quite reliable for electronic states in 
high Landau Icvels^^'^^. The situation should be simi- 
lar for graphene, particularly if one can show that LLM 
is small for states in high LLs, as we will indeed find 
self-consistently below. We thus adopt the Hartree-Fock 
approximation for the states we study. 

More specifically, our Hartree-Fock approach to the 
Dirac equation description of uniform and stripe quan- 
tum Hall phases in graphene is adapted from a procedure 
developed for electrons in a standard 2DEC^^; in what 
follows we briefly outline the method, and highlight the 
(largely technical) differences. The Hamiltonian for the 
system in a magnetic field is 



(2) 



1234 



where the numbers denote composite indices for the 
different quantum numbers specifying the states [e.g., 
1 = (ni, ATi, si, ti) = (LL index, guiding center coor- 
dinate, spin, pseudospin (valley) index)]. 



hvp 

T 



ei = sgn(ni)^^V2K| - si5*^i3-B (3) 



is the LL spectrum plus the Zeeman energy, and 

^^1234 = jWi234Ss,sJt,tJs,sJt,t, Hi^)'"-" (4) 



i=l 



are matrix elements for the Coulomb interaction. W1234 
is related to standard matrix elements^^ 



i^V(q)(ni,Xi|e"^-nn4,X4) (5) 



S 



x{n2,X2\e-"i-^\m,X3), 
with V{q) — 27re^/(7 and 
{ni,Xi\e"i-'-\7n,Xi) 



exp[i-fe(-'^i + X4^)]F„,^n4iQ)Sxi,Xi+qyi^, 



where 



X exp 



T n\ —714 

4 / V 2 
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for n4 < ni, where L'^{x) is the generahzed Laguerre 
polynomial. Note that F„^,„i(q) = [F„i^„J-q)]*. 
In terms of V, W takes the form 

+sgn(nin4)V|„j|_ij„2|^|„3|j„^|_i 
+sgn(n2n3)Vj„j|j„^l_ij„3|_i^l„^l 
+Sgn(nin2n3ri4)y|„^l_i^l„2|_ij„3|_ij„^l_i]. 

Note that the guiding center coordinates {X) have been 
suppressed in the subscripts in Eqs. ([5]) and jB]). The 
density matrix operators are defined as 

I 



This relation may be inverted to obtain the expectation 
value of an arbitrary single-particle operator in terms of 
density operator expectation values, 



P 

X exp[^ipx{X + X')]Sx,x'-pye.^- 



For states with discrete translational symmetry, the sum 
over p is restricted to reciprocal lattice vectors {G}. The 
interaction part of the HF Hamiltonian Hhf may now 
be written as 



"2, "3 G S2,t2 Si,ti 

where 

UH{n2,n,;G) = ^ E E ^^'("i' ^3, 714; G) (p^J^^J^i (-G)) , (10) 

rii .714 si,ti 

Ux{si,S2,ti,t2,n2,n3;G) = — ^ ^sC^i' "2, na, 714; G) (-G)) , (11) 

ni ,714 

with 

4 

== (_,)l«il+l"2|,|n3|+|n4| J](V2)*-«[i?(|ni|, \n,\, \n2\, |n3|; G) + sgn(nin4)i/(|ni | - 1, \n,\ - 1, \n2\, \ns\;G) 
1=1 

+sgn(n2n3)_ff(|ni|, |n4|, \n2\ - 1, |n3| - 1;G) + Sgn(nin2n3n4)i/(|ni | - 1, |n4| - 1, \n2\ - 1, |n3| - 1;G)], 
-'^g("-i,n2,n3,n4;G) 

4 

= (_j)l"il+l"2ljl«3|+|n.| J|(y2)^-.o[X(|ni|, |n3|, |n2|, |n4|;G) + sgn(nin3)X(|ni | - 1, |n3| - 1, \n2\, \ni\;G) 

i=l 

+sgn(n2n4)X(|ni|, jngl, l^al - 1, |n4| - 1; G) + sgn(nin2n3n4)X(|ni | - 1, {n^l - 1, \n2\ - 1, |n4| ~ 1;G)], 

where 

ff(ni,n2,n3,n4;G) = ^^y(G)F„„„,(G)F„3,„,(-G), (12) 

X(ni,n2,n3,n4;G) = — ^ ^(9)-^ni,n2(9)-Pn3,™4(-9) exp(-iq x Gf). (13) 

q 

The single-particle Green's function is defined by 

G;:?*-,- = - (rc„;,,,(r)cl,^,,,,,(0)) , (14) 
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and its Fourier-transform by 

G-*,,,(G,r) = '^-^Y.G-^Jl,AX,X ^ Gye,T)eM-^G..X +]^G^Gye). (15) 

Within the HF approximation, the equation of motion (EOM) for GJJf*,j, (G, w™) is given by 

(zc^,„ + A*/?i)G:^vt'(G,^m) -\ Ag{s,t,n,si,h,n3;G,G')G'^f:^,'MG',uj„,) - Snn'Sss'Su'Sc^o, (16) 

ti,Si,ra3,G' 

where 

Ag{s,t,n,si,ti,n3;G,G') = e'l^Sn^nStttSstsScG (17) 

+ [UH{n, ng; G' - G)6t,t6s,s - Ux{si, s, h,t, n, n^; G' - G)]e^^><«'^'/^ 



Because LLM could be important, we retain several "ac- 
tive" LLs (with LL indices between niower and riupper) 
around the chemical potential (see Fig. i.e., we solve 
the EOM explicitly for the Green's function matrix al- 
lowing off-diagonal elements in the LL index for values 
n satisfying niower < n < Uupper- However, it would be 
incorrect to completely neglect the filled LLs below the 
active LLs (ric <n< niower)- These levels can enter the 
calculations through Uh and Ux ■ However if sufficiently 
below the chemical potential, we expect LL mixing to be 
negligible for these states. We thus treat these levels as 
"inactive" , and fix their density matrix elements to be 
{p^fl,^,{G)) = Snn'Sss'Stt'Sco- (For self-consistency, we 
verify numerically that LL mixing for the lowest active 
level is very small, justifying the dividing point between 
active and inactive levels.) With this form the inactive 
levels do not contribute to Uh due to the (1 — 5g,o) in the 
Hartree term; i.e., it is precisely cancelled by an interac- 
tion with a uniform neutralizing background. However, 
these levels do contribute a non-vanishing exchange en- 
ergy Ux, 

Ur'\si,S2,h,t2,n2,n3;G) = 



Sx('^2, "3)^5132^*1*2^,0, 



where inact stands for "inactive" and 
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Sx(?^2,?^3)= ^ Xg{ni,n2,n3,ni]0). (18) 



FIG. 2: Division of LLs into active and inactive LLs. In 
this example, niower = 1 and riupper = 5. The up- 
per and lower cutoffs of the active LLs are determined by 
the self-consistency conditions 1 — {pn''°"""^sti^)) ^ 1 ^^'^ 
^Pn"pperst(0)^ ^ 1- Chemical potential /i for the case of 
filling = 14 is indicated as an example. In general sev- 
eral active levels are retained both above and below n for all 
calculations reported here. 



We can rewrite Eq. ([TT]) as 

Ag{s,t,n,si,ti,n3;G,G') 

= [^n^nsn- J^^X {n, n3)]St^tSsisSG' G 

+ [U^^'{n,n3;G' -G)St,tSs,s 

-Un.su s,h,t,n,n3;G' - G)]e»«x«'^'/2, 

where the superscripts act means now the summations in 
Uh and Ux are restricted to the active LLs. 



The calculations involve solving Eq. (|T6|) to obtain the 
Green's function, from which we obtain the density op- 
erator matrix elements. Finally, the Hartree-Fock energy 
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is given by 

Ehf 

S 



E E{[<.'^"-3-^sx("2,n3)](p-:?(o)) 

"2. "3 s,t 

G 

-Y.Ur{3i,S2,h,h,n2,ns;G){pZlltAG))]} 



Si ,t 

s 



EE 



where the last hne is a constant for given B and nia 



TABLE I: The diagonal density matrix elements (pnstiO)), 
indicating the occupation of the spin and pseudospin split LL 
with quantum numbers (nst). The occupations are very close 
to a step function, indicating that LLM is weak. The devia- 
tion from a step function decreases as B decreases, indicating 
that LLM becomes less important. Here g is set to 1. 



n 


St 












He = 1872 


= 12000 


= 24000 






(B = 2or) 


{B = 3.12r) 


{B = 1.56r) 



5 (iJ;),(ifr),(T4) 0.0000697757 0.000059254 0.0000559402 

5 (Tfr) 0.000520925 0.000425574 0.000396389 

4 (i4),(ifr),(Tl|) 0.00078437 0.000651856 0.000611152 

4 (Tfr) 0.499997 0.50002 0.500026 

3 (iJ;),(ifr),(TJ|) 0.999181 0.999318 0.999361 

3 (Tfr) 0.999567 0.999627 0.999647 



III. RESULTS OF THE CONTINUUM LIMIT 
MODEL 

Table U details some typical results for the occupa- 
tions of the various Landau levels near the Fermi en- 
ergy. In this example v = 14.5; i.e., the LL with 
{nst) = (4 Iff) is half-filled. Note that in this notation 
we denote the valley index as a pseudospin, with two val- 
ues and IJ- denoting the K and K' valleys respectively. 
We present results for different coupling constants in the 
range 0.5 ^ g < 1, consistent with previous estimates 
of its appropriate valueiii^. Our qualitative results are 
very similar for different values of g, even for (unphysi- 
cal) values well above 1. We can see that the occupations 
immediately become very small above the half-filled LL, 
and very close to 1 below it, indicating that LLM is in- 
deed small. This small level of mixing, in spite of the 
small non-interacting energy gap between LLs where the 
Fermi energy is located, may be understood as being a 
consequence of the large exchange enhancement of the 
gap due to the filled LLs. Furthermore, for smaller B, 
deviations of the occupations from a step function de- 
creases (albeit just slightly), which means for decreasing 
field LLM becomes even less important. This unintu- 
itive result occurs because of the large sea of negative 
energy LL states. With smaller field the degeneracy of 
each of these decreases, and so that more inert LLs need 
to be filled to obtain the correct density of electrons [see 
Eq. ID)]. In units of e^/ei, the exchange interaction in- 
creases with increasing \nc\, and the LLs effectively be- 
come slightly more separated. 

Fig. [3] illustrates the LLM for two integer fillings where 
the system is in a uniform liquid state. Here we measure 
the LLM by the quantity 

where the sum is over active LLs only {niower — —5 and 
n-upper = 5). We again see that LLM is small and de- 
creases as B decreases for fixed filling factor. 
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FIG. 3: LLM for integer fillings. Note ric oc 1/B, indicating 
the LLM decreases with decreasing field. 



Previous studies of crystal and stripe states in 
graphene in which a single^^ d^d^ or small number^^ of 
Landau levels is retained find that such states can be sta- 
ble in the presence of a magnetic field. Our study sug- 
gests that inclusion of the large number of LLs intrinsic 
to graphene not only does not change such results, but 
even increases their validity in weak fields. The result 
of this is that, within a zero-temperature mean-field de- 
scription, one expects that in a very clean system many 
different states will persist down to the origin in a phase 



diagram plotted in the n vs. B plane [see Fig. 1(b) )]. The 
state is determined only by the filling factor. This is in 
sharp contrast to the situation for conventional 2DEG's, 
where LLM always destabilizes such states as the origin 
is approached. 

One seeming paradox associated with this behavior is 
how the system approaches the uniform state which is be- 
lieved, at least within a mean-field approach, to occupy 
the origin of the n vs. B phase diagram for graphene. 
The answer lies in noting that since LL mixing is negligi- 
ble, the relevant length scale for the charge-ordered states 
of a partially filled LL is the magnetic length, which di- 
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FIG. 4: Wavelength and amplitude for stripe with u — 14.5, 
niower ~ 1 and Uupper ~ 5, as functions of the cutoff ric (and 
alternatively of the field B). The insets are density profiles for 
different g's, with the arrows indicating the quantities plotted. 



verges as i? 0. Fig. [5] illustrates the consequence of 
this for stripe states. One sees that the wavelength and 
amplitude of the density modulation are basically con- 
stants when measured in appropriate units {£ and l/27r£^, 
respectively). Thus, these quantities should, up to loga- 
rithmic corrections, follow simple scaling relations, 



wavelength cc £ cc , 
amplitude oc XjP' oc B. 



(19) 
(20) 



As B decreases, the stripes, and we believe CDWs in 
general, "evaporate" , and thus approach the expected 
uniform density state at the origin. 
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FIG. 5: Geometry of the tight binding problem. The unit 
cell corresponds to the area between the two armchair chains, 
including one of the bounding chains. In this example the 
unit cell contains 8 sites and 4 plaquettes. In our calculations, 
there are Uy unit cells in the y direction, but only one in the 
X direction. We apply periodic boundary conditions to both 
X and y directions, so that site 8 is connected to site 1 by the 
bond to its right, etc.. 



IV. HUBBARD MODEL 

The continuum limit forces one to adopt a cutoff in the 
occupied states, which in the previous section was accom- 
plished by adopting an appropriate choice of the mini- 
mum occupied LL index, Uc- Since the increase of Uc with 
decreasing field tends to suppress LLM, one may wish to 
consider whether a more physical cutoff scheme would 
give similar results. Towards this end we re-examine this 
question within a tight-binding Hubbard model. As in 
the continuum case, we look for states of this system 
within the Hartree-Fock approximation. For a simple 
on-site interaction U, the HF Hamiltonian for spin up 
electrons is 



E 

<ij> 



' all flit' 



(21) 



where < ij > indicates nearest neighbors. For spin down 
electrons the Hamiltonian is analogous, with "f and | in- 
terchanged. 

We choose the unit cell to be the area between two 
adjacent armchair chains (see Fig.[S]). We apply periodic 
boundary conditions to both x and y directions and study 
stripes oriented along the y directions. We can Fourier 
transform along the y direction, then we only need to 
define the phases of tij along one armchair chain (e.g., 
the chain 12345678 in Fig. [5]), i.e., — > ii.i+i, where 
i = 1,2, . . . jTix labels the sites in the unit cell {rix = 8 
for the example in Fig. [5] and in general can be any integer 
multiple of 4). One possible choice is 



argi 







(-l)('-i)/2(i- l)7ra 



if i is even 
if i is odd 



where a — with $m being the total number 

of flux quanta in the unit cell. (See Table |TT1) Since 
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TABLE II: Nonzero argti,i+i 


for hopping from site i 


to i + 1 


on the chain 12345678 in Fig 


El 




i 3 


5 


7 




Ana 


— 67ra 



<i>m must be an integer, the magnetic fields for compu- 
tationally tractable system sizes are actually very large 
compared to experiments. Nevertheless, we can deduce 
the qualitative behavior from these calculations. 

The coupling constant in this model is g = Ui/ta, 
where t « 2.7 eV is the hopping amplitude of the tight- 
binding approximation. The situation that g is field- 
independent does not arise naturally here; we introduce 
it by adjusting U with field according to the relation 

C/cx VScx J—. 

V Tlx 

For real Coulomb interactions, the effective HF po- 
tential includes a short-range exchange potential and a 
long-range Hartree potential, both proportional to ^/B. 
Stripe and bubble states result from the competition of 
these?^^-'^^. Because of the highly local nature of the inter- 
action in the Hubbard model, neither this scaling nor the 
effective long-range part of the interaction emerge: one 
only finds a local repulsion between electrons of different 
spins. Thus, the charge-ordered bubble and stripe states 
are not eigenstates of Eq. (PTjl : one generically finds uni- 
form density states. To obtain the former states, one 
needs to include longer-range interactions in the Hamil- 
tonian. Obtaining full solutions of the HF approximation 
in this situation is possible but challenging, and is unnec- 
essary for our more modest goal of testing the effect of 
using a real lattice rather than an energy cutoff. Thus, 
rather than fully solving for states of a system with long- 
range interactions, we include a slowly varying external 
potential which models the effect of the long range (non- 
contact) part of the potential. For simplicity we take this 
to have the form 

Ai?(t) = A'^cos(2Tix,/L^)a\^a,ii, 

i 

where A must scale with field in the same way as U , and 
Lx is the length of the unit cell along the x direction. 

Our goal is to study how the density of a CDW state 
varies if the field is allowed to change, keeping the effec- 
tive g fixed. In order to make a fair comparison between 
states at different field strengths, we also fix the ratio 
^ oc ^/rix^m so that the width of the stripes and their 
spacing relative to the unit cell size does not change. This 
restricts the number of systems we can examine. How- 
ever, the data we do get are in excellent agreement with 
the continuum model, i.e., the stripe amplitude (defined 
as the difference in maximum and minimum densities) is 
roughly proportional to the magnetic field (see Table IIIII 
and Fig. [S]). Note that in these calculations the ampli- 
tude decreases slightly faster than linearly with the field. 



0.015 



•S 0.010- 
^ 0.005 - 



0.000 ^ — ' — ' — • — ' — ' — ' — • — ' — ' — ' — • — ' — " — ' 
0.000 0.002 0.004 0.006 0.008 

— cc B 

FIG. 6: Stripe amplitude for Hubbard model calculation 
for a fixed ratio of unit cell width to magnetic length, as 
a function of ^m/rix, which is proportional to the field. 
U/t = 76.2^$™/n:,, A/U = Q.l. A straight fine through 
the data points extrapolates rather accurately through the 
origin. 

TABLE III: Change in stripe amplitude, defined as the differ- 
ence between the maximum and minimum densities within a 
unit cell, when the magnetic field is changed by varying "I>m 
and Ux- For first row of data, A/U — 0.2; for second row, 
A/U = 0.15. For all data, U/t = 76.2y^¥^Jn^. 



($m, ria:) 




B'/B 


amplitude' /amplitude 


(1 , 600) 


(2 , 300) 


4 


4.1941 


(1 , 800) 


(2 , 400) 


4 


4.1509 



consistent with a decreasing role for Landau level mix- 
ing. This effect is larger for larger values of A/U, as 
illustrated for example in Table Hill 

In Fig. [S] we illustrate the stripe amplitude for states 
generated for three values of ($m, J^k), corresponding to 
three different magnetic fields, but with the ratios of the 
unit cell sizes and magnetic length the same, and with a 
relatively small value of A/U (0.1). In this case one may 
fit a straight line through these points, and find that it 
extrapolates to the origin rather accurately. This is con- 
sistent with the stripe amplitude continuously vanishing 
in the _B — > limit, as was found in the continuum ap- 
proach. 



V. SUMMARY 

We have examined the stability of liquid and charge- 
ordered states for graphene (focusing on stripes as a 
paradigm for the latter) in the quantum Hall regime 
against the effects of Landau level mixing. Because the 
coupling constant g is independent of field, we find the 
LLM does not increase with decreasing field, and that, 
counterintuitively, it decreases, albeit by a small amount. 
This latter effect is due to a large exchange enhancement 
of the LL gaps from the filled negative energy LLs, which 
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increase in number as the field decreases. Within mean- 
field theory, this implies that clean and cold graphene 
at small fields and densities should support many differ- 
ent phases, determined solely by the filling factor. This 
contrasts with the conventional 2DEG, where a Wigner 
crystal state is believed to reside throughout this regime. 
In graphene, the liquid phase thought to exist in the ab- 
sence of doping is reached in the B limit at fixed 
filling factor by an "evaporation" of the CDW, in which 
the amplitude vanishes linearly with B. 
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